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Abstract. In the largest solar energetic particle (SEP) events, acceleration occurs at shock 
waves driven out from the Sun by coronal mass ejections (CMEs). Peak particle intensities are 
a strong function of CMP- speed although the intensities, spectra, and angular distributions of 
particles escaping the shock are highly modified by scattering on Alfv&i waves produced by the 
streaming particles themselves. Element abundances vary in complex ways because ions with 
different values of Q/A resonate with different parts of the wave spectrum, which varies with 
space and time. Just recently, we have begun to model these systematic variations theoretically 
and to explore other consequences of proton-generated waves. 

INTRODUCTION 

Before the discovery of CMEs, it was believed that all SEPs were accelerated in so- 
lar flares. This left serious questions of how the particles could cross from one mag- 
netic field line to another to yield high intensities over ~18(f of solar longitude and why 
events lasted lor days when the flares themselves were over in a few hours at most. 
With the observation of CMEs and improving observations of SEPs, it became clear 
that the largest “gradual’' proton events were associated with CME driven shock waves 
(11). At the highest energies in the largest events, measurements of the mean charge 
state of Fe found that Qr. - 14 at energies out to 600 MeV/amu (46). This showed 
that these ions came from ambient coronal temperatures and not from flare-heated ma- 
terial. Furthermore, these ions would be rapidly stripped at coronal densities (40) and 
were probably accelerated beyond 2 solar radii. A distance of >5 solar radii was inde- 
pendently derived by Kahler ( 10) for peak acceleration of 1 - 21 GeV protons observed 
in ground-level events. 

Meanwhile, SEPs that do come from impulsive flares are the smaller “ ? He-rich” 
events (27, 32) where Q ic ~ 20 ( 17). The unique 1000-foki enhancement of ’He/^He in 
these events is believed to come from resonant wave-particle interactions in the flare 
plasma (e.g. 39). 

Most physicists practicing today learned the origin of SEP events from their re- 
search advisers in simpler times when the “solar flare myth” (7) prevailed. Perhaps be- 
cause research advisors are often as revered as gods, the idea seems to have lodged as 
a religious belief that is difficult to shake, and shock acceleration is still viewed as her- 
esy in some quarters. However, those willing to tempt fate will find that a consistent 



picture of distinct “gradual" and “impulsive" events emerges that spans a wide range of 
observations from photons to ionization states to abundances and spectra from 20 keV 
to 20 GeV. 

A complete comparison of the properties of gradual and impulsive SEP events and 
their sources is beyond the scope of this paper (see reviews 27, 29. 30. 32, 14. 15. 9, 
7). Recent multi-spacecraft studies of the spatial distribution of SEPs around the CME 
and shock have been reported by Rcames, Barbier. & Ng (33) and by Reames. Kahler, 
& Ng (34). These include the phenomenon of spatially and temporally invariant spectra 
produced when particles are quasi-trapped in expanding magnetic bottles formed by 
CME loops or behind the shock (see also 32). 

This paper seeks only to highlight certain features of SEP events. It is not a review. 
For a more comprehensive review of both gradual and impulsive SEP events and their 
relationship to other heliospheric acceleration sources, see the recent review by Rcames 
(32). In the present paper, we consider two recent updates on the properties of CME- 
associated SEP events, energy-dependent variations in ionization states and correla- 
tions between SEP intensities and CME speeds. We then discuss numerous conse- 
quences of wave-particle interactions rear 
shocks and finally consider the high-energy 
“knee" of the spectrum where shock accel- 
eration efficiencies plummet. 

IONIZATION STATES 

One of the distinctive signatures of shock 
acceleration has been the ionization states of 
the elements from 0.3 to 600 MeV/amu. es- 
pecially Fe, that are similar to those of the 
solar wind. Not only do these ionization 
states represent ambient, unheated coronal 
plasma, but, in addition, the ions must be ac- 
celerated at high altitudes above the corona 
where densities are too low for stripping to 
occur. 

Figure 1 shows new measurements of 
ionization states in the event of 1997 No- 
vember 6 (18. 20). These indicate strong en- 
ergy dependence in just the region where the 
ion speed begins to exceed the thermal elec- 
tron speed. If shock acceleration begins low 
enough in the corona that the density is ~10 9 
cm' 3 , ions might come to an equilibrium 
charge state appropriate to their speed rela- 
tive to the plasma electrons. A simple exten- 
sion of the common expression for equilib- 
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FIGURE 1. Mean ionization state of ele- 
ments in the 1997 November 6 SEP event 
(18, 20) are compared with a simple theory 
of stripping during shock acceleration in 
the upper corona (38). 



rium charge states (38) is compared with the data in Figure 1 . The measured Fe charge 
at 40 MeV/amu falls below' the equilibrium curve, perhaps indicating either that Fe has 
not traversed enough material to attain equilibrium or that acceleration continued after 
the shock left the corona, or both. 

Within the present accuracy of the measurements, the energy dependence of the 
charge states in this large SEP event Is consistent with shock acceleration in the low- 
corona. At low energies, when the ion speed is less than the thermal electron speed, 
ionization states still define a plasma temperature that is consistent with an ambient 
coronal temperature of ~2 MK, 

One should be aware of recent papers, such as that of Barghouly and Mewaldt (1). 
that completely neglect ion speed and determine ionization states using the plasma 
temperature alone, even for ion energies up to ~100 MeV/amu! This paper leaves the 
erroneous impression that shocks are incapable of producing highly ionized ions. 

CORRELATIONS OF PEAK INTENSITIES AND CME SPEEDS 

Another strong indication of the origin of particles in the large SEP events is the 
correlation between peak proton intensities and maximum CME speed that was origi- 
nally observed by Kahlcr et al. (11). We have recently reexamined this correlation for 
two different samples of events (35). The first list of events was used in a recent ISTP 
workshop with CME speeds determined by the SOHO/LASCO coronagraph, and par- 
ticle intensities measured on the Wind and IMP-8 spacecraft. The second is an older list 
of events using the SOI. WIND coronagraph at Earth while the Helios spacecraft ob- 
served protons off the solar limbs. Both lists were prepared to look at correlations of 
CMEs and interplanetary shock waves, without regard for the presence or absence of 
energetic particles. Peak proton intensities, or upper limits, were determined at ~ 2 and 
~2() MeV for these events and their correlation with CME speed Is shown in Figure 2. 

The correlation coefficients (r=0.62 and 0.72) listed in Figure 2 are quite high 
compared with corresponding correlations (r~0.3) of particle intensities with flare 
properties such as X-ray or y-ray intensities. In fact, the correlation with CME speed is 
surprisingly high, considering that many factors are not considered, such as the longi- 
tude of the event, the CME speed on the observers actual field line, the speed of the 
material upstream of the shock, and the angle between the magnetic field and the shock 
normal and its time evolution. 

Because of the steep dependence of particle intensities on CME speed, only the fast- 
est 1-2% of all CMEs drive fast enough shock waves to accelerate any particles at all. 
Most CMEs are emitted with speeds that barely exceed that of the solar wind. This 
speed correlation strongly suggests that it is the shock ahead of the CME where accel- 
eration occurs, rather than the reconnection region behind it, for example. The exis- 
tence of many large, slow CMEs (-98% of them) with a complete absence of energetic 
panicles further supports this conclusion. 
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FIGURE 2. Correlations of peak proton intensities at -2 and ~20 MeV with CME speed are 
shown for 2 different event lists. Correlation coefficients are relatively high and LS fit lines 
show a steep dependence of the intensities on speed (35). 

PROTON GENERATED WAVES 

It is a fundamental fact of plasma physics that protons streaming along magnetic 
fields can generate or amplify Alfven waves (41). These waves modify the transport of 
particles that come behind. Different particle species at the same velocity resonate with 
different regions of the wave spectrum because of their different values of charge-to- 
mass ratio. Q/A. Hence, element abundances can vary in space and time in response to 
the changing wave spectrum. The role of self-generated waves in acceleration at inter- 
planetary shocks was first considered in detail by Lee (13). These waves produce what 
is referred to nowadays as a “cosmic-ray modified” shock. 

Ironically, waves that can scatter particles are notoriously difficult to resolve from 
the background turbulence by direct measurements. Fisk (6) pointed out long ago that 
particle transport is grossly inconsistent with observed magnetic power spectra. Pre- 
sumably this occurs because much of the variation lies in a direction orthogonal to the 
magnetic field and is not seen by the particles traveling along the field (2). However, 
the enhanced spectra of sell-generated wave near shocks, where streaming is largest, 
have been found to be consistent with the expectations of Lee (13) theory (12. 42. 43. 
49). Of course, the newly-generated wave energy at a fast shock is less than a few per- 
cent of the energy in SEPs. which is smaller than the magnetic energy. 



Streaming-Limited Intensities 

Figure 3(a) shows superposed intensity-time profiles of several events observed on 
Helios (26). Early in these events proton intensities do not exceed -200 (errr sr s 
MeV)' 1 . although later they may rise 1 or 2 orders of magnitude at the time of shock 
passage. There are no events near 1 AU with significantly higher intensities of MeV 
protons early in the event. Ng and Reames (23) showed theoretically that these intensi- 
ties were limited by self-generated waves produced at a shock near the Sun when 
streaming protons produced sufficient wave intensities to throttle the streaming. In suf- 
ficiently large events, higher-energy particles are also streaming-limited, but at lower 
intensities as shown in Figure 3(b). Here again, peak intensities occur at the time of 
shock passage, even at 100-500 MeV. Recently, Reames and Ng (37) undertook a sta- 
tistical study of the streaming limit, its energy dependence and its radial gradient. 
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FIGURE 3. Panel (a) shows superposed intensity-time profiles of 3-6 MeV protons in several 
events with streaming-limited intensities early in the events. Panel (b) shows similar limits as 
a function of energy in the large 1989 October 19 event. Note that intensities often peak at the 
time of shock passage (32). 

Abundance Variations 

Measurements of energetic particles at a few MeV/amu have been used to determine 
element abundances in the solar corona for many years (22. 3. 28, 31, 32). When the 
element abundances are averaged over many events, are divided by the corresponding 
photospheric abundances, and are plotted as a function of the first ionization potential 
(FIP), one obtains the distinctive plot that is believed to describe the ion-neutral frac- 
tionation during formation of the solar corona. 


In dealing with event-to-event variations, workers in this field have historically con- 
sidered all large events. Fe-rich to Fe-poor. either ranking them by their Fe/O ratios 
(22) or averaging over all the events (3. 28. 31). Breneman and Stone (3) found that 
element abundances in an event had a rough power-law dependence on Q/A of the ele- 
ments, both for Fe-rich and Fe-poor events. They were able to rank ~20 elements in 
this way. and to determine coronal abundances for each. Had Breneman and Stone 
used modem photospheric abundances for their comparisons, they would have found 
that the abundances, averaged over all of their events, could be directly interpreted as 
“coronaf abundances, without correction. 

However, two elements that failed to fit the Q/A dependence were simply omitted 
by Breneman and Stone, namely H and He. Surely, one must be wary of a model that 
fails to include the dominant elements of the plasma. The elements C through Zn cor- 
relate well with Fe; based on 49 large SEP events, for example. Si/O vs. Fe/O gives a 
correlation coefficient of 0.90 (28). However, He/H vs. Fe/O gives a coefficient of 
only 0.31 (28). Plots of abundances vs. Q/A show disorganized behavior when H and 
He are included (32). One can hardly blame this behavior on an uncertainty about the 
ionization state of H. Our discipline has harbored this secret about H and He for many 
years. Yet, things would only get worse before they got better. 

For many years there had been some evidence for time dependence of abundances 
during an event, but we were unprepared for the smooth systematic variations we be- 
gan to see with the large-geometry in- 
struments on the Wind spacecraft. 

These variations are seen in Figure 4, 
where the abundances in the lower panel 
are normalized to coronal abundances 
(see 47, 48). 

A physical basis for understanding 
the abundance variations was provided 
by applying the same theory of self- 
generated waves that clarified the 
streaming-limited intensities and spectra 
Ng et at. (24, 25) modified the earlier 
calculations of particle transport through 
self-generated waves (23) by injecting 
power-law energy spectra with coronal 
abundances at the position of the shock 
as it moved out toward Earth. Wave 
generation and pitch-angle scattering of 
the ions was followed in full detail in 
space and time. Abundance variations 
from this theory (24, 25) are shown in 
Figure 5 for conditions similar to those 
of the 1998 April 20 west -limb event. 
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FIGURE 4. Selected intensities and abundance 
ratios, relative to coronal, are shown for the 1998 
April 20 event. Abundances respond to the vary- 
ing wave spectra ahead of the shock. 



The behavior of the variations in 
Figure 5 may be understood qualita- 
tively as follows. Karly in the event. 

Fe is scattered slightly less than O of 
the same velocity, even though the 
scattering mean free paths. of 
both species are ~1 AU in the pre- 
sumed Kolmogorov, i t . back- 
ground wave spectrum. Hence, Fe 
arrives first and Fe/O begins high but 
falls as more C) arrives. As wave 
growth increases with time. 0 is 
preferentially trapped near the shock 
and Fe/O at 1 AU rises. To obtain 
the final decease in Fe/O and the 
greater decrease at higher energy, 
we must decrease the shock com- 
pression ratio with time, as might he 
quite appropriate for an event near 
the west limb of the Sun. For a 
shock that remains stronger. Fe/O 
would remain high and its value would increase with energy as is seen in other events. 

The He/H ratio in Figure 5 is especially interesting. Its behavior is almost a mirror 
image of that of Fe/O. However, early in the event, before the first protons to arrive 
have had time to grow waves, one might expect He/H to behave like Fe/O. Both in- 
volve the ratio of high-rigidity/low-rigidity species at the same velocity; both would 
respond similarly to a Kolmogorov wave spectrum. The secret lies in the fact that He 
resonates with waves created by protons of twice its velocity. These waves are created 
prior to the arrival of the He and scatter the 2. 1 MeV/amu He while the 2. 1 MeV pro- 
tons have yet to generate resonant waves. Thus, in large events , with substantial wave 
growth, He/H can rise initially; in smaller events, the ambient power-law wave spectra 
will produce an initial decline in both Fe/O and He/H. The energy region of the initial 
rise in He/H (and other species) can expand to higher energies in larger, more intense 
events (24, 25). The physical origin of the initial rise or fall in He/H has been a mystery 
to observers for the last 20 years (50). 

With the above discussion, we also begin to understand the poor correlation be- 
tween He/H and Fe/O on an event-averaged basis. Furthermore, at its present level 
the theory only considers differential transport of ions, not differential acceleration. If 
Fe/O is enhanced in one place, it must be reduced in another; the spatial integral of all 
species should yield coronal abundances. Thus, it becomes clearer why averaging 
abundances over many events, Fe-rich and Fe-poor, recaptures coronal abundances and 
produces the familiar plot of averaged abundance ratios as a function of FTP. 



0.2 


0 10 20 30 40 

Tim© (hrs) 

FIGURE 5. Abundance variations calculated from 
the theory of Ng ei ai (24) show the same qualitative 
behavior as the observations shown in Figure 4. 
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FIGURE 6. Energy thresholds of various instruments are 
superposed on Fe/O data from Tylka et al. (44). The gradual 
events in 1989 have Q Fe ~ 14 up to 600 MeV/amu (46); Q Ft 
for the 1997 Nov. 6 was discussed and is shown in Figure 1 . 


Based on recent abun- 
dance measurements at sig- 
nificantly higher energy on 
ACE/SIS, Cohen et al. (4) 
have suggested that while 
Fe-poor events may come 
from shock acceleration. 

Fe-rich events come from 
Hares, or worse, some other 
unknown mechanism. 1 
would suggest that the true 
problem here is that the 
measurements were made at 
an energy where large 
abundance variations are 
known to exist. Tylka et al. 

(44) previously found en- 
hancements of Fe extending to much higher energies, even in the same events of Sep- 
tember and October of 1989 where they found Qpt-14 from 200-600 MeV/amu (46). 
These Fe/O ratios are shown in Figure 6 together with the regions of energy coverage 
of the ACE/SIS, Wind, and ISEE-3 instruments. Furthermore, Mazur (19, 21) studied 
the energy dependence of abundances in 10 large SEP events. He found large abun- 
dance variations at high energies while, at low energies, the event-to-event abundance 
variations become smaller and the mean abundances approached those of the corona 
and solar wind. All of the abundance measurements of Reames (28, 31,32) were made 
using ISEE-3 data At a few MeV/arnu the resolution of impulsive and gradual events 
is clear (36), above 15 MeV/amu, it is not. If we examine the 4 events studied by 
Cohen et al. (4). the abundances seen on the Wind spacecraft at -2 MeV/amu. in the 
same events, are less dramatic and are consistent with the historic values for gradual 
events. The high-energy abundances are also consistent with observations during the 
last two solar cycles; unfortunately, Cohen et al. neither discuss nor reference this ear- 
lier work. Ionization-state measurements confirm that the large Fe-rich SEP events are 
from the same CME-driven shock sources as are the Fe-poor ones. Therefore, it is ap- 
propriate to average over all large SEP events to obtain abundances. The fact that this 
averaging directly produces coronal abundances is no accident. 

Angular Distributions 


Another prediction of the Ng (24, 25) model is that wave generation in sufficiently 
large SEP events rapidly reduces scattering mean free paths, not only at the shock, but 
also out as far as 1 AU. Figure 7 compares intensity-time profiles and magnetic- 
sectored intensities of H and He in “large” and “small” events seen by the Wind space- 
craft. The sectored data are produced by binning individual ions identified onboard 
into 16 sectors spanning 360° relative to the direction of the magnetic field as the 
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FIGURE 7. Intensities and angular distributions relative to the magnetic field arc for compared 
for “small” and “large" SEP events. Wave generation in the large events rapidly reduces stream- 
ing along the field and isotropizes the particles. 


spacecraft spins with spin axis normal to the ecliptic plane. The magnetic azimuth is 
updated every 3 seconds onboard the spacecraft and high-speed onboard processing 
can accommodate all particles that enter the EPACT/LEMT telescope. Figure 7 shows 
initial particle flow centered on 180° to the field direction early in both events. How- 
ever. this flow persists for more than a day in the small 1998 June 16 event, but the 
particles rapidly isotropize within hours in the 1998 April 20 event. Peak proton inten- 
sities differ by a factor of ~100 in the two events at ~2 MeV; evolution of the scattering 
mean free paths for the April event simulation are show by Ng et al. (24). 

This wave -generated scattering in large SEP events explains another historic con- 
troversy. Early measurements, which were limited to large events by the instrument 
sensitivity, found scattering mean free paths. A<(). 1 AU, white later observations of 
small impulsive-flare events often found A>1 AU (see review 32). Could the inter- 
planetary medium distinguish the particles from different events and scatter them differ- 
ently? Ironically, the answer is yes lor events that occur at different times. Protons 
from large events rapidly modify the interplanetary medium; they are no longer “test 
panicles” that probe pre-existing turbulence. 


The “Knee” of the SEP Spectrum 

After a sufficiently long time, the spectrum of shock-accelerated particles will begin 
to approach an equilibrium power-law form extending to arbitrarily high energies with 
an index that depends only on the shock compression ratio (8). For real shocks, how- 



ever, this equilibrium is not attained at 
high energies where the scattering 
mean free path and the acceleration 
time may be long. 

Lovell ei al. ( 16) showed a proton 
energy spectrum extending from 1 
MeV to 10 GeV for the large SEP 
event of 1989 September 29 as shown 
in Figure 8. This event is a well- 
studied shock event from a source just 
beyond the west limb of the sun; it is an 
event for which Qp c -14 was observed 
in the 200-600 MeV/amu region (46). 
The spectral region from 1 to 10 GeV. 
shown as shaded in Figure 8. was ob- 
tained by fitting ground-level neutron- 
monitor measurements from many sta- 
tions. The fitted curves for H and He 
in the figure come from the shock the- 
ory of Ellison and Ramaty (5). which 
has the form of a power law times an 
exponential in energy/nucleon. The ex- 
ponential rollover in this theory, which 
scales as Q/A. was originally used to ex- 
plain differences in electroa proton and 
He spectra in large SEP events. 

In smaller events, the spectral “knee” 
occurs at much lower energies. Ele- 
ment abundances near the knee in the 
1 998 April 20 event are shown in Figure 
9 (45). In this event, the e-fokiing con- 
stant, Eo. in the exponential rollover 
scales like Q/A times energy/nucleon. 

Comparison of Figures 8 and 9 raises 
an obvious question. Why does the pro- 
ton spectral knee occur at -20 MeV in 
one event and at -1 GeV in another? 
This large difference can have little to do 
with ambient scattering conditions for 
high-energy protons prior to the event. 
It may occur because oik event has a 
faster shock, higher proton intensities, 
and much greater wave growth that can 



Figure 8. Energy spectra in the 1989 September 
30 event as measured on IMP 8, GOES, and by 
ground-level neutron monitors (shaded region) 
(16). Fitted curves are from shock theory (5). 
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FIGURE 9. Energy spectra of various ion spe- 
cies in the 1998 April 20 event (45) observed by 
IMP 8, WIND, and ACE. 



persist to much higher energies. The position of the knee Is a consequence of the 
strength of the shock and the intensity of proton-generated waves. 

CONCLUSIONS 

New measurements of particle intensities, ionization states, spectra and abundances, 
together with new observation of CMEs, continue to support the premise that accelera- 
tion in large SEP events occurs at CME-driven shock waves. Of course, observations 
in some energy regions are too highly variable to be helpful in resolving these gradual 
events from impulsive- flare related events. 

Many of the features we observe in the large SEP events are a coasequence of pro- 
ton-generated waves near shocks. These include the following: 

1) Streaming-limited intensities early in large events (23. 24, 26. 37. 32). 

2) flattened low-energy spectra in large events (24, 32. 47). 

3) Systematic time variations in abundances (31, 32. 24. 25. 47, 48). 

4) Abundance variations that average to coronal (HP-dependent) values (32). 

5) He/H uncorrelated with Fe/O; breakdown of the power law of abundances vs. 
Q/A resulting from non-Komolgorov wave spectra (24. 25. 32. 47). 

6) Initial rise in He/H in large events (24. 25, 47. 32). 

7) Rapid onset of isotropy in large event s, even at 1 AU. 

8) Large variations in the energy of the spectral “knee” (32, 45). 

Amazingly, all these results arise from wave energy that is less than a few percent of 
the energy in SEPs. This list has come from the close collaboration of observations 
and theory that is as uncommon as it is highly rewarding. Not only has the new theory 
been able to explain historic questions that have lain dormant up to 20 years. It has 
also predicted new behavior that was subsequently observed. 

CME-driven shock waves provide a richly diverse laboratory in which to study the 
physics of particle acceleration. Shock acceleration occurs widely at sites throughout 
the heliosphere. However, only in SEP events can we examine that process for hun- 
dreds of events under such a wide range of initial conditions. We are only beginning to 
organize the diversity we see. 
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